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Parameter-free and perturbatively stable leading order (LO) and next-to-leading order (NLO) parton densities 
for real and virtual photons are presented. 



1. Introduction 

The partonic content of real and virtual pho- 
tons can be measured, for example, in electron 
positron scattering. The measured e + e~ cross 
section can be obtained by a convolution of a flux 
of target photons B with the cross section of deep 
inelastic electron photon scattering 



da(ee — > eX) 



/ 7 (pa)/ e * da(ej(P 2 ) -> eX) 
+0(P 2 /Q 2 ) . (1) 



The virtuality of the target photon, P 2 , has to 
be much smaller than the virtuality Q 2 of the 
probe photon such that the neglected terms of 
order 0(P 2 /Q 2 ) are small. 

In the following, we are interested in the pho- 
ton structure function \x, Q 2 ) which enters 
the deep inelastic electron photon cross section in 
exactly the same way as is well known from deep 
inelastic electron proton scattering 



dcr(e 7 (P 2 ) -> eX) 
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F%^ P ' is a convolution of the massless parton 



densities of the real or virtual photon with the re- 
spective massless Wilson coefficients, in order to 
obtain a factorization scheme independent, i. e., 
physically meaningful result 



Fl {P \x,Q 2 ) = 

u,d,s 

a 

+ 2n 



MP 2 



+ 



— e 2 C 
27r e <? G 7 



(3) 



This also holds for the direct photon coefficient 
C^, i. e., C 7 has to be the massless direct photon 
coefficient (irrespective of P 2 ) since it is related to 
the massless photon to quark and photon to gluon 
splitting functions Pg 7 and P$ in the evolution 
equations. 

As usual, we work in NLO in the DIS 7 scheme 
in which the destabilizing log(l — x) of C 7 , as 
calculated in the MS scheme in (||), is absorbed 
into the quark distributions 
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Again, it should be emphasized that we use the 
real photon coefficient irrespective of P 2 . 

2. Boundary Conditions 

The parton distributions are then obtained by 
QCD evolution of appropriate boundary condi- 
tions at a low scale Qq rj 0.3 GeV 2 , which will 
be presented in the following. The exact LO 
and NLO values of the universal (i. e. hadron- 
independent) input scale Qo are fixed by the 
experimentally well constrained radiative parton 
densities of the proton (3j . 

2.1. Real Photon 

The boundary condition for the real photon ||] 
is given by a vector meson dominance (VMD) 
ansatz where (at the low scale Qo) the physical 
photon is assumed to be a coherent superposition 
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Figure 1. Comparison of our radiatively generated 
LO and NL0(DIS 7 ) predictions for F% (x, Q 2 ), based 
on the valence-like parameter-free VMD input in eq. 
(^), with the data of ref. |J. For comparison the 
GRV 7 results are shown as well. In both cases, the 
charm contribution has been added, in the relevant 
kinematic region W > 2m c , according to fixed order 
perturbation theory. 



of vector mesons which have the same quantum 
numbers as the photon 

P(x, Ql) = f had (^ Ql) = <*G 2 f r°(x, Ql) (5) 

with G 2 ud = (g p ± g^) 2 and G 2 g = G 2 S = g 2 p + gl 
{g 2 p = 0.50, gl = 0.043). This optimal coher- 
ence maximally enhances the up quark which is 
favoured by the experimental data. 

Since parton distributions of vector mesons 
(p, u,...) are unknown we furthermore assume 
that these are similar to pionic parton distribu- 
tions. Thus, for given pionic parton distributions 
our model has no free parameter. 



2.2. Pion 

Since only the pionic valence density is exper- 
imentally rather well known, we utilize a con- 
stituent quark model to relate the pionic light sea 
and gluon to the much better known parton dis- 
tributions of the proton || . In this model the pro- 
ton and pion are described by scale independent 
constituent distributions (U p , If 7 * ,...) convoluted 
with their universal partonic content f c [Q 
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where / = v, q, g. 

Since in Mellin-n-space the convolution is a 
simple product, the ratio / f p is independent 
of the flavour / and one can easily find boundary 
conditions for the pionic gluon and sea which only 
depend on the rather well determined valence dis- 
tribution of thepion and the parton distributions 
of the proton 
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^<f, - q «{n,Ql) = -- q 

2.3. Virtual Photon 

In this subsection boundary conditions for the 
much more speculative virtual photon are pro- 
posed. 

First recall that the photonic parton distribu- 
tions can be written as a sum of a pointlike and 
a hadronic part 



7(P 2 ) 



(x,( 



x,Q 2 ).(8) 



The pointlikc part is a particular solution of the 
inhomogenous evolution equations of the pho- 
ton and vanishes by definition at the input scale 
P 2 = max(P 2 ,Q2): f^ p2 \x,P 2 ) = 0. While 
the pointlike solution is perturbatively calcula- 
ble, the homogenous (hadronic) solution requires 
a boundary condition. At P 2 we assume the 
hadronic part of the virtual photon to be given 
by the hadronic part of the real photon sup- 
pressed by a rho-meson propagator r](P 2 ) = (1 + 
P 2 /m 2 p )- 2 : 



'had 
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This boundary condition of course smoothly ex- 
trapolates to the real photon case. 

Note that, the employed dipole suppression 
factor is somewhat speculative and can be re- 
garded as the simplest choice of modelling the 
P 2 -suppression. 
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Figure 2. Detailed small-a; (as well as large-x) be- 
havior and predictions of our radiatively generated 
u 1 — u 1 and g 1 distributions in LO and NLO(DIS 7 ) 
at fixed values of Q 2 . The dashed-dotted curves 
show the hadronic NLO contribution /^ ad to / 7 . The 
valence-like inputs at Q 2 = Q 2 = Mlo.nlo ||j ac- 
cording to eq. (^), are shown by the lowest curves 
referring to fi 2 . For comparison we show the steeper 
NLO GRV 7 expectations as well. The results have 
been multiplied by the number indicated in brackets. 



3. Numerical Results 

3.1. Real Photon 

In Fig. [I] we compare our LO and NLO predic- 
tions for F2 , given by the full and dashed lines, 
with experimental data [0. Also shown are the 
GRV 7 H predictions which are steeper in the 
small- a; region (due to the steeper sea, which has 
been generated from a vanishing input). In all 
cases the charm contribution has been added ac- 
cording to a fixed order calculation (see, e.g., eqs. 
(15) and (16) in @). 

As one can see, we achieve a good description 
of the data, particularly in the region of small Q 2 
where the GRV 7 and also the SaS ID || expec- 
tations fall somewhat below the data. Further- 



more, our LO and NLO results are close together 
demonstrating an excellent perturbative stability. 

Figure || shows the x-dependence of the up 
quark and the gluon for three values of Q 2 . The 
solid and dashed lines are our NLO and LO re- 
sults, respectively. The dotted lines are the NLO 
GRV 7 distributions. 

Comparing the solid lines with the dashed- 
dotted lines, representing the hadronic contri- 
bution to the solid lines, one can see that for 
x > 0.1 the pointlike part dominates (especially 
for the quark distribution) , whereas at very small 
x the hadronic component is dominant implying a 
small-a; behaviour similar to the pion or proton. 

Also shown are our valence-like inputs which 
are (vanishingly) small for x < ICC 2 . Thus, the 
small-a; increase at higher scales is purely due to 
the QCD evolution. 
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Figure 3. Comparison of our LO and NLO pre- 
dictions for xg 1 at Q 2 = (Mf) 2 ) = 75 GeV 2 with 
HERA(Hl) measurements juj. The GRV 7 and SaS 
expectations are taken from refs. Q and M, respec- 
tively. 



In Fig. |3j we compare the a:-dependence of our 
gluon distribution with recent HI measurements 
pC|| . The full data points are dijet events and the 
open points are high pt track events. Our LO 
and NLO results are very similar to the gluon of 
SaS ID, shown by the dotted line. The flatter 
small-a; behaviour compared to GRV 7 is partly 
induced by the flatter gluon distribution in the 
GRV-98 proton ||. (Another reason is that the 
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GRV 7 gluon is enhanced by a factor k = 1.6(2.0) 
in NLO (LO) at Q 2 , in eq. (2) in §].) 
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Figure 4. LO and NLO (DIS 7 ) predictions for the 
up-quark and gluon distributions of a virtual photon 
7(P 2 ) at Q 2 = 10 and 100 GeV 2 . For comparison 
the results for the real photon (P 2 — 0) are shown as 
well. For P 2 = 0.5 and 1.0 GeV 2 the NLO 'hadronic' 
contribution in (H) is also shown separately. The GRS 
expectations are taken from ref. |llj . The results have 
been multiplied by the numbers indicated in brackets. 



3.2. Virtual Photon 

Figure [| is the virtual photon analogue of Fig. 
|^. Again the cc-dependence of the up quark and 
the gluon is shown for two values of Q 2 and sev- 
eral photon virtualities P 2 . The solid and dashed 
curves are our NLO and LO distributions, respec- 
tively, whereas the dotted curves are the NLO 
predictions of Gliick, Reya and Stratmann [TT| , 
denoted by GRS. For Q 2 = 100 GeV 2 we also 
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Figure 5. Predictions for the LO effective parton 

density xf l{p2) (x,Q 2 ) at the scale Q 2 = (p^ = 

85 GeV 2 and at two fixed values of x. The HI data 
Jl2| have been extracted from DIS ep dijet produc- 
tion. The solid curves refer to our predictions in the 

theoretically legitimate region P 2 <C Q 2 = (Pj? 
whereas the dashed curves extend into the kinematic 
region of larger P 2 approaching Q 2 where the concept 
of parton distributions of virtual photons is not valid 
anymore (see text). 



show the LO GRS up quark distributions (wide 
dotted), which are very similar to our LO results. 

On can see that the pointlike part becomes 
increasingly important with increasing P 2 due 
to the dipole suppression of the hadronic com- 
ponent [r?(P 2 (GeV 2 ) = 0.5,1.0,5.0,10.0) = 
0.3,0.14,0.01,0.003]. At P 2 = 5,10 GeV 2 the 
hadronic part is nearly vanishing even at small 
x. Furthermore, note the perturbative stabil- 
ity of the parton distributions guaranteeing the 
required stability of the physical structure func- 
tions. 

Finally, in Fig. |^ we compare our LO parton 
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distributions with an effective parton density 



q—u,d,s 
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(10) 



which has been extracted from HI dijet data [|f 2| . 
Q = jfy* and we show results for two values of 
x in dependence of the photon virtuality. The 
transition from the solid to the dashed line in- 
dicates that for P 2 approaching Q 2 the concept 
of renormalization group resummed parton dis- 
tributions of virtual photons is not expected to 
hold any more, because the resummed logarithms 
log Q 2 /P 2 arising in the partonic subprocesses are 
not much larger than the non-logarithmic terms. 

In the relevant kinematic region P 2 « Q 2 the 
agreement with the data is reasonably good. This 
is in contrast to a simple dipole suppression of the 
real photon //(P 2 )/ 7 ^^ ) (x, Q 2 ), as illustrated 
by the dotted curves in Fig. g, which fails to de- 
scribe the data. 
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